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Abstract 

A method is presented for determining the initial conditions of clas- 
sical orbits from the quantum spectra of the diamagnetic hydrogen 
atom. Each classical trajectory which is closed at the nucleus pro- 
duces a sinusoidal fluctuation in the photoabsorption spectrum. The 
amplitude of each orbit's contribution appears in the Fourier transform 
of a spectrum computed at constant scaled energy. For a given initial 
state, closed-orbit theory gives the dependence of this recurrence am- 
plitude on the initial angle of an orbit. By comparing the recurrence 
amplitudes for different initial states, the initial conditions of closed 
classical orbits are determined from quantum spectra. 
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1 Introduction 

According to the Correspondence Principle, the solutions to Schrodinger's 
equation for a given system should contain the classical motion of the same 
system as a subset (the "classical limit"). In practice, however, gaining 
information about the classical dynamics of a system from the quantum 
behavior has proved challenging, particularly for chaotic systems. In some 
cases, the distribution of eigenvalues can be used to discern whether the 
classical motion is regular or chaotic [HE]- Spectral periodicities are related 
to periodic orbits [3J H] and can be used to determine such classical quantities 
as the period, action, and stability of the periodic orbits [HE]. This paper 
shows that spectral periodicities of the diamagnetic hydrogen atom can be 
used to determine the initial conditions of the periodic orbits which pass 
through the origin (closed orbits). 
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Semiclassical quantization techniques reverse the causal role between 
quantum and classical behavior by using the classical solutions to construct 
approximate quantum solutions. EBKM tori quantization [3] gives approx- 
imate eigenvalues of conservative integrable systems. Periodic-orbit theory 
predicts the level density from a sum over periodic orbits and is applicable 
in regimes of integrable and non- integrable classical motion [3] . The closely 
related closed-orbit theory [6] predicts photoabsorption spectra from a sum 
over the orbits which are closed at the nucleus. 

Periodic-orbit theory can also be used to glean classical information from 
quantum level densities [5]. The main result of periodic-orbit theory [3] is 
an expression for the fluctuating part of the level density, 

oo 

9c(E) = E E A nk e i <- nS >- a <**' 2 \ (1) 

k n=l 

(atomic units). The index k distinguishes the primitive periodic orbits: the 
shortest period orbits for a given set of initial conditions. S/% is the action of 
a periodic orbit. Each primitive orbit retraces itself, leading to new orbits 
with action nSk, where n is an integer. Hence, every repetition of a periodic 
orbit is another periodic orbit. The quantity A nk is related to the stability 
of an orbit, and a n k is the orbit's Maslov index. 



2 Closed Orbit Theory in Diamagnetic Hydrogen 

The Hamiltonian of diamagnetic hydrogen is 

H = V - - - + -B 2 p 2 , (2) 
2 r 8 H W 

where the magnetic field B is taken to be along the z axis. This Hamiltonian 
can be scaled so that for spectra at constant scaled energy, e = EB' 2 / 3 , the 
fluctuating part of the level density is [S] 

oo 

9cH = E E A nk e^ nw ^- a ^' 2 \ (3) 

k n=l 

where w = B~ 1 / 3 , and S = S/(2ttw) is the scaled action. At fixed e, 
the classical dynamics has no dependence on w. Consequently, the Fourier 
transform of an energy level spectrum computed at fixed e as a function of 
w gives peaks which lie at the scaled action of periodic orbits and whose 
heights are related to the stability of the orbits [7j. 
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Closed-orbit theory [6] is similar to periodic-orbit theory, except that 
closed-orbit theory is applicable only to atomic and molecular spectra and 
yields the oscillator strength density from a specified initial state whereas 
periodic-orbit theory yields the density of states. Only orbits that begin and 
end at the nucleus are important in closed-orbit theory. Physically, these are 
associated with the outgoing waves that are generated when a tightly bound 
electron is excited to a high-lying state. For diamagnetic hydrogen, every 
orbit which is closed at the nucleus is also a periodic orbit whose period is 
equal to either the closure time or twice the closure time. 

According to closed-orbit theory, the average oscillator strength density 
at constant e is given by a smooth background plus an oscillatory sum of 
the form [6] 

oo 

f(w) =^2^2 D l nk sm(2>irnwS k - cf) nk ). (4) 

k n=l 

4>nk is a phase that depends on the Maslov index and other details of the or- 
bit. D l nk is the recurrence amplitude of a closed orbit for a given initial state 
(labeled i). It contains information about the stability of the orbit, its initial 
and final directions, and the matrix element of the dipole operator between 
the initial state and a zero-energy Coulomb wave. The Fourier transform 
of an oscillator strength spectrum computed at fixed e as a function of w 
is called a recurrence spectrum, because it gives peaks which correspond to 
the scaled action of closed orbits and whose heights correspond to D l nk . 

3 Method for Determining Initial Angles 

The dependence of the recurrence amplitude D l nk on the initial and final 
angles of an orbit can be used to determine these angles from computed 
photoabsorption spectra. We can write the recurrence amplitude as 

D i nk = F i (9 i J f )G nk , (5) 

where Oi and f are the initial and final angles the orbit makes with the z axis 
at the origin. (For a specified scaled energy, the initial conditions of a closed 
orbit are completely determined by its initial angle.) The dependence on 
the initial and final angles of the orbit and on the initial state is completely 
contained in the function F t (9i,9f). G n k is related to the stability of the 
orbit [6\. Several factors simplify the discussion for the scaled energy (e = 
—0.7) considered here [8]. The dynamics is near integrable in this regime 
and all closed orbits have either Of = Oi or Of = tt — 0j. Furthermore, 
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Fig. 1: Predicted ratios of recurrence amplitudes from different initial states. 



F l (6i,9f) = F l (6i,iT — 9 f). Consequently, the angular dependence can be 
considered a function of a single angle, 6. 

The form of F{6) for a given initial state can be obtained from Du and 
Delos [6j. For m = final states, 

F 3s (9) = 0.625 sin 6>[Pi(cos#)] 2 (6) 
F 2p (9) = 0.5797sm9[P 2 (cos6)+P (cos9)/4\ 2 (7) 
F 3d {9) = 0.8505 sin 9[P 3 (cos9) + 2P 1 (cos9) /Q] 2 , (8) 

where the Pi{cos9) are Legendre polynomials. (For orbits parallel to the 
field, there is no sin (9 term.) The angular functions are normalized [9] such 
that 

r F i (9)d9 = l. (9) 
Jo 

Since all the dependence on the initial state is contained in F l (9), 

Dk = m. (10) 

Consequently, the ratio F l (9)/F°(9) for a particular closed orbit can be de- 
termined by dividing the recurrence amplitudes for the corresponding peak 
in the computed recurrence spectrum. 

Fig. 1 shows F 2 p(9)/F 3s (9) and F 3d {9)/F 3s (9). These ratios are sym- 
metric about 90°. A given ratio corresponds to either two or four possible 
initial angles, but since a closed orbit of diamagnetic hydrogen with initial 
angle 9 has a symmetric partner of the same shape with initial angle tt — 9, 
only the range from < 9 < tt/2 is of interest. This leaves either one or 
two possible initial angles, depending on the value of the ratio, but this 
ambiguity can be resolved by looking at both ratios. 
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Fig. 2: Recurrence amplitude for odd-parity, m = final states at e = —0.7 excited 
from the 3s initial state. (The recurrence amplitude is the square root of the power 
spectrum.) 



4 Recurrence Spectra 

We now turn to the recurrence spectra. A number of methods are available 
for computing spectra of diamagnetic hydrogen [2j [10|. Illj . Diagonalizing 
the Hamiltonian matrix in a spherical hydrogenic basis is adequate for e = 
—0.7 [12] . The oscillator strength spectrum, f(w), is multiplied by w 3 to 
remove the global variation of oscillator strength over the range of fields 
(31 < w < 133) and normalized so that 

/"VJ max 

/ w 3 f(w)dw = l. (11) 

This produces narrower peaks in the recurrence spectrum and allows for 
comparison with recurrence spectra of different initial states. (In princi- 
pal, one can obtain arbitrarily high resolution in the recurrence spectrum 
by computing spectra to very large w. This allows identification of orbits 
arbitrarily close in action.) 

The m = odd-parity recurrence spectrum for the 3s initial state is 
shown in Fig. 2. (van der Veldt et al. have analyzed a similar experimental 
recurrence spectrum for the 2s initial state in helium [13].) To interpret 
the spectrum, it is helpful to recall the basic structure of closed orbits in 
diamagnetic hydrogen. At large negative scaled energies, there are three 
primitive short-period orbits: one moves on the p axis, and two move on 
the ±z axis pp. Most longer-period primitive closed orbits are created by 
bifurcations of these orbits and their repetitions. 

The orbits on the p axis and those which bifurcate from them are called 
rotators. Peaks corresponding to these orbits will be labeled i?* , where n 
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Fig. 3 Top: Recurrence amplitude for even-parity, m = final states at e = —0.7 
excited from the 2p initial state. Bottom: Recurrence amplitude for odd-parity, 
m = final states at e = —0.7 excited from the 3s initial state. This recurrence 
spectrum is multiplied by —1 for comparison. 



denotes the repetition of the orbit on the p axis from which the orbit bifur- 
cated, and b distinguishes between different orbits which bifurcated from the 
same parent. The orbits on the z axis and those which bifurcated from them 
are called vibrators. Recurrences corresponding to these orbits are labeled 
where n and b have the same meaning as for the rotators. For S < 4 
the recurrence spectrum is dominated by repetitions of the perpendicular 
and parallel orbits (labeled R n and V n , respectively) [H] . For S > 4 other 
orbits are present which were born by bifurcations of these orbits, and the 
recurrence spectrum becomes more complicated. 

5 Determination of Initial Angles 

This method for determining initial angles can be tested by considering 
recurrence spectra for 9.9 < S < 12.4, which are shown in Fig. 3 for the 2p 
and 3s initial states. This range of action has a number of orbits with non- 
trivial (not 0° or 90°) initial angles and is sufficiently low that the density 
of orbits allows most recurrences to be resolved. 

The method for determining the initial angle of an orbit is demonstrated 
by considering the peak at S = 10.0661. The amplitudes for the 2p and 3s 
initial states are D 2p = 0.0361 and D 3s = 0.0294. Their ratio corresponds 
to F 2p (9)/F 3s (6) = 1.228. Two initial angles give this ratio: 9 = 28.51° and 
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Table 1: Predicted initial angles from two different amplitude ratios are compared 
with actual classical initial angle, cl . All angles are in degrees. 
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9 = 73.47°. This ambiguity can be resolved by considering the amplitude for 
the 3d initial state. D 3d = 0.0181 which corresponds to F 3d (6) / F 3s (6) = 
0.615. Once again, two initial angles give this ratio: 9 = 29.07° and 9 = 
53.90°. Consequently, the initial angle is determined to be near 9 = 28.51°. 
This is in good agreement with the actual initial angle, 9 cl = 29.10°. 

This method can be used to determine the initial angle of the closed 
orbits corresponding to many of the peaks in Fig. 3. Table Q] shows that this 
method gives agreement with the initial angles obtained classically. Some 
of the closed orbits corresponding to peaks in Fig. 3 are shown in Fig. 4. 

The initial angles with the largest errors are for orbits where one or more 
of the recurrence amplitudes is small or where the ratio F l (9) / F 3 (9) has a 
small slope. If the recurrence amplitudes themselves are small, there is a 
larger relative error in the computed amplitudes due to noise in the Fourier 
transform. There is also a larger relative error in the semiclassical approx- 
imations for small recurrences, because the semiclassical approximations 
used to derive Eq. [5] neglect the contribution of a closed orbit's neighbors 
to the angular dependence of its recurrence strength [14]. The error in the 
predicted initial angle is large if the ratio of angular functions has a small 
slope because a small error in the computed recurrence amplitudes or in the 
semiclassical angular functions leads to a large error in angle. In principle, 
the errors in the semiclassical approximations become smaller as Ti — * 0. 
This limit can be approached practically in diamagnetic hydrogen by look- 
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Fig. 4: Some closed orbits corresponding to peaks in Fig. 3 and Table [TJ (a) V^; 
(b) (c) Rf 5 ; (d) R\ 5 . 

ing at smaller magnetic fields at a given scaled energy. The difficulty of 
small slopes in the angular ratios can be remedied by considering a range of 
initial states so that, for every angle, at least one of the angular ratios has 
a large slope. A weighted average over many initial states should produce 
more accurate initial angles than consideration of a few initial states. 

6 Conclusion 

In summary, a method is presented for determining the initial conditions 
of classical orbits from the quantum spectra of the diamagnetic hydrogen 
atom. However, this quantum mechanical method is an inefficient choice 
for searching for closed classical orbits. More efficient classical methods are 
available for finding closed and periodic orbits by directly integrating Hamil- 
ton's equations [15] . This development of a quantum mechanical method of 
finding closed orbits is a practical step toward understanding how the solu- 
tions to Schrodinger's equation yield the solutions to Hamilton's equations 
in the classical limit. It has been shown to be accurate in the near- integr able 
regime of diamagnetic hydrogen. A simple generalization should be appli- 
cable in more strongly chaotic regimes. 
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